Ferroelectricity in ultra-thin perovskite films 
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We report studies of ferroelectricity in ultra-thin perovskite films with realistic electrodes. The 
results reveal stable ferroelectric states in thin films less than 10 A thick with polarization normal to 
the surface. Under short-circuit boundary conditions, the screening effect of realistic electrodes and 
the influence of real metal/oxide interfaces on thin film polarization are investigated. Our studies 
indicate that metallic screening from the electrodes is affected by the difference in work functions 
at oxide surfaces. We demonstrate this effect in ferroelectric PbTi03 and BaTi03 films. 



The effect of size on thin-film ferroelectricity has been 
known for a long time, but has not been completely un- 
derstood. Initial experiments and mean field calculations 
based on the Landau theory suggested that below a crit- 
ical correlation volume^ of electrical dipoles between 10- 
100 ran 3 , ferroelectricity vanishes due to intrinsic size 
effectsi2i& For thin films with the polar axis perpendic- 
ular to the surface, incomplete compensation of surface 
charges creates a depolarizing field that has been shown 
to further reduce the polarization stability4£ 

Recently, however, monodomain ferroelectric phases 
have been observed in very thin films, below ten unit- 
cells thick. 6 - 7 - 8 Furthermore, Fong et at showed that fer- 
roelectric phases can be stable down to ~12 A (three 
unit cells) in PbTiC>3 films by forming 180° stripe do- 
mains, suggesting that no fundamental thickness limit is 
imposed by the intrinsic size effect in thin films^ This 
idea has been corroborated by ab initio calculations car- 
ried out on perovskite films, which tell that no critical 
thickness exists for polarization parallel to the surface^ 
and that polarization perpendicular to the surface can 
exist in films three unit-cells thick if the depolarization 
field is artificially removed piLiS 

On the other hand, it has been found that the depo- 
larization field plays a dominant role in reducing polar- 
ization normal to the surface and depressing ferroelectric 
transition temperatures in thin films. In a continuum 
model by Batra et alA which includes the depolariza- 
tion effect, a critical thickness of 100 A for perovskite 
films was analytically derived, assuming a Thomas-Fermi 
screening length of 1 A for the metal electrodes. Also, 
a recent first-principles calculation revealed that BaTiC>3 
films with SrRuC>3 electrodes lose ferroelectricity below 
^24 A (6 unit cells), 13 thus suggesting that a minimum 
thickness limit exists for useful ferroelectric films. While 
indeed a minimum film thickness must be influenced by 
the polarization of the ferroelectrics and the screening 
length of the electrodes^ it is not yet clear whether 
the depolarizing field can ever be completely removed 
by realistic electrodes on ultrathin films, nor how mon- 
odomain thin film ferroelectricity is affected by the choice 
of electrodes and by the interactions at the metal/oxide 
interface. 

In this paper, we provide answers to these questions. 



In particular, we investigate how the critical thickness 
varies in different systems with realistic electrodes. We 
find that even though ferroelectricity is lost in BaTiOa 
thin films, polarization close to the bulk value can be 
stabilized in PbTiC>3 thin films with thickness less than 
10 A; the behavior in BaTiC>3 is therefore not universal 
to all ferroelectric thin films. Furthermore, we show that 
inequivalent ferroelectric/electrode interfaces can assist 
in stabilizing ferroelectricity in thin films. 

We apply density functional theory (DFT) calcula- 
tions^ to ultrathin ferroelectric capacitors that are con- 
structed of PbTiC>3 and BaTiC>3 films sandwiched be- 
tween two conducting electrodes. Two electrodes com- 
monly used in ferroelectric devices are studied: plat- 
inum (Pt) and the metallic oxide SrRu03. Short-circuit 
boundary conditions are imposed by the periodic bound- 
ary conditions and electrodes of sufficient thickness. 

We consider AO (PbO or BaO) and Ti02 ferroelectric 
terminations, and we examine different ferroelectric- 
electrode interfaces, focusing on the lowest energy 
one for each termination. Pt is most stable with 
the first layer situated above the oxygen atoms on 
the Ti02 terminated surface, and above A and O 
atoms on the AO-terminated surfaced The period- 
ically repeated supercells can thus be described by 
the general formula^ Pt 4 /AO-(Ti0 2 -AO) m /Pt 5 and 
Pt 4 /Ti02-(AO-Ti0 2 ) m /Pt5 with Pt electrodes and 
(SrO-Ru0 2 )2/AO-(Ti02-AO) m /(Ru02-SrO)2-Ru0 2 
and (Ru0 2 -SrO) 2 /Ti0 2 -(AO-Ti0 2 ) m /(SrO-Ru0 2 ) 2 - 
SrO with SrRuOs electrodes (periodic boundary 
conditions means that the ferroelectric slabs are sepa- 
rated by nine layers of electrode). Figure ^ shows the 
structures for two representative systems at m = 2. 
Although ferroelectric instabilities can be present in the 
directions parallel and perpendicular to the film, here 
we focus only on the latter situation, as our primary 
interest is the depolarization effect. Therefore, for all 
the capacitors considered, the in-plane atomic positions 
are kept fixed at the ideal perovskite positions and the 
in-plane lattice constants are set equal to the experi- 
mental value for the corresponding bulk ferroelectric 
perovskite, i.e. a=3.935 A for PbTi0 3 , and a=3.991 A 
for BaTi0 3 i£ 

We start by determining the fully relaxed structures 
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FIG. 1: (a) Ti0 2 -terminated Pt/PbTi0 3 /Pt structure, (b) 
PbO-terminated SrRu0 3 /PbTi0 3 /SrRu0 3 structure. The 
ferroelectric displacement patterns are shown by arrows. The 
rumpling of each layer (displacement of cations relative to 
anions) is marked, in A. 



TABLE I: Polarization (P) and tetragonality ratio c/a for 
PbTi0 3 films with Pt and SrRu0 3 electrodes. The bulk val- 
ues are given in the last line. 



termination electrode 
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Pt 
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0.89 
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1.00 
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bulk 


0.75 


1.060 



of the supercells. The atomic positions and strains are 
fully relaxed along the surface normal direction. Fig- 
ure Wi. a ) shows the relaxed supercell structure of the 
Ti02-terminated PbTiC>3 film with Pt electrodes and 
m=2. The displacements from the centrosymmetric 
structures are calculated. The relative displacement be- 
tween the cations and anions for each PbTi03 layer (see 
Figure HJa)) clearly demonstrates that the structure is in 
a ferroelectric state. As TableUshows, PbTiC>3 films with 
Pt electrodes are ferroelectric for all thicknesses down to 
m=l. For both terminations, the polarization values in 
these systems are slightly larger than the corresponding 
bulk valued 

These results suggest that bulk ferroelectric polariza- 
tion can be stabilized in thin films below 10 A with real- 



istic electrodes. As Figure^illustrates, an upward point- 
ing polarization leads to an enhancement of the displace- 
ments at the bottom interface, and a reduction at the top 
interface, relative to the interior layers. This observation 
is in agreement with previous DFT calculations with ex- 
ternal fields^ In going from m=2 to m=4, both P and 
c/a decrease towards the bulk values as the surface-to- 
volume ratio decreases and the surface effect is averaged 
over more layers. To further check whether there exists 
a thickness limit below which ferroelectricity disappears 
in this system, we examine the structure at m — 1. We 
find that a single unit-cell has a stable polarization of 
0.89C/m 2 , despite the stoichiometrically different envi- 
ronment as compared to bulk PbTi03>& We therefore find 
strong evidence for an absence of a critical thickness for 
ferroelectricity in PbTiC>3 films with Pt electrodes. In 
contrast to the PbTiC>3 films, none of the BaTiC>3 struc- 
tures were found to be ferroelectric for m=2 or m=4, in 
agreement with the results in Ref. 1 13 , 

To make contact with earlier analytic thcory^^ we 
examine the electric field in the Pt/PbTiOs/Pt capaci- 
tor by plotting the macroscopic-averaged^ electrostatic 
potential, as shown in Figure [3 Also shown in Figure [21 
is the potential in the free-standing PbTiOa film in which 
a bulk ferroelectric displacement perpendicular to the 
surface is imposedi^l The slope of the potential in the 
PbTiC>3 slab, which is the depolarizing field, is signifi- 
cantly smaller in the capacitor (0.009 V/A), than in the 
free-standing slab (0.324 V/A). Our calculations show 
that in the absence of electrodes, the latter field brings 
the system back to the paraelectric structure. The can- 
cellation of a substantial fraction (97%) of the depolariz- 
ing field is due to metallic screening from the grounded 
electrodes that compensates the polarization charge. 
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FIG. 2: Macroscopic-averaged electrostatic potential 
along the film normal direction of the Ti0 2 -terminated 
Pt/PbTiOa/Pt capacitor (solid curve) and a free-standing 
PbTiOa slab with a fixed bulk ferroelectric displace- 
ment (dashed curve). The potential is constant inside the Pt 
electrodes. 

In addition, the screening is accompanied by the for- 
mation of unequal local dipoles at the two interfaces due 
to different chemical bonding, as shown in Figure [3J On 
the top interface, the Pt and Ti atoms lose charge, which 
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is redistributed between the atoms, forming a Pt-Ti alloy. 
On the bottom interface, where the Pt-0 distance is the 
shortest, the Pt and O atoms lose charge, the Pt from 
the d z 2 orbitals, and the O from the p z orbitals^ while 
the Pt d xz and d yz orbitals gain charge. Similar behav- 
ior is observed at the ferroelectric/electrode interfaces of 
the AO-terminated capacitors. This inequivalent charge 
arrangement at the two ferroelectric/metal interfaces is 
consistent with the different interface polarizations we 
noted earlier. 
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FIG. 3: Contour plot of the induced charge density for the 
Ti02-terminated Pt/PbTiOa/Pt capacitor structure at the 
(a) top interface (b) bottom interface. (The Pt layers farther 
from the interfaces are not shown.) Electron loss is given by 
solid lines and electron gain by dotted lines. The induced 
charge density is found by subtracting the charge densities 
of free-standing Pt and PbTiOs slabs from the total charge 
density in the capacitor. 



With SrRu03 electrodes, the heterostructures also 
adopt a ferroelectric state, as shown in Figure^b) . How- 
ever, unlike with Pt electrodes, the polarization with 
SrRu03 electrodes is only half the bulk value, indicating 
that the surface charges are only partially compensated. 
The c/a ratio for the SrRuOs/PbTiOs/SrRuOs capacitor 
falls to 1.045 as the film thickness decreases to two unit 
cells. The corresponding electrostatic potential shows a 
depolarizing field of 0.07 V/A across the PbTiOs slab, 
which is much smaller than in the free-standing PbTiOs 
film and explains why the polarization is stabilized at 
a significant level. Nevertheless, the stronger field com- 
pared to that in the Pt/PbTiOs/Pt structures also con- 
firms the weaker screening effect of SrRu03 compared to 
Pt and other transition metals^ 

We note further that allowing the internal coordinates 
of the SrRu03 electrodes, especially those in the bound- 
ary layers, to relax together with the PbTiOs ions is cru- 
cial. When the SrRuOa ions are fixed in their ideal po- 
sitions, we find only paraelectric structures for either of 
the PbTiOs terminations. This result indicates that the 
interaction between the PbTiOs and the SrRu03 elec- 
trodes plays a crucial role in stabilizing ferroelectricity: 
the screening effect alone cannot account for the signifi- 
cant reduction of the depolarizing field. 

The results presented so far have suggested that fer- 
roelectric polarization normal to the surfaces can be sta- 
bilized in thin films less than 10 A thick. However they 



do not explain why the ferroelectric instability can be 
retained in PbTiOs but not in BaTiOs films when the 
same electrodes are applied. To address this question, 
we turn now to look at the depolarizing field and screen- 
ing effects in these systems. The depolarizing field was 
previously addressed in phenomenological studies r>"i ,2 $ 
and is shown to scale with the spontaneous polarization 
of the ferroelectrics, the screening length of the electrodes 
and the inverse ferroelectric film thickness. Hence the de- 
polarizing field would be expected to be negligible only 
in very thick films. 

In our calculations, inspection of the electrostatic po- 
tential in the free-standing film with a fixed bulk ferro- 
electric displacement shows that the potential drop across 
the ferroelectric slab (Ai) is different from the potential 
difference between the two asymptotic vacuum potentials 
(A2), as illustrated in Figure^] This difference arises be- 
cause the two surfaces have different work functions, as a 
result of the polarization orientation, parallel to the top 
surface normal and anti-parallel to the bottom one. Be- 
cause Pauli repulsion keeps metal electrons out of the fer- 
roelectric, the potential drop that is "seen" and screened 
by the electrodes is A2, not Ai. An electrostatic anal- 
ysis 2 ^ of the potential and the electric field in the slabs 
shows that this difference must be taken into account 
when modeling the screening in all realistic systems. 
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FIG. 4: The electrostatic potential of a free-standing film 
with a fixed bulk ferroelectric displacement (lower curve); the 
self-consistent potential from solving the Poisson equation for 
Thomas- Fermi screening charges (upper curve). Panel (a) is 
PbTi0 3 film and panel (b) is BaTi0 3 film. 

We use a simple Thomas-Fermi model to highlight the 
effect of metallic screening on a ferroelectric slab, with- 
out the complex structural and electronic relaxations of 
DFT calculations. We treat the DFT macroaveraged 
electrostatic potential of the free-standing film as an ini- 
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tial potential and find the metallic screening charge us- 
ing Thomas-Fermi theory. We then find the total self- 
consistent potential by solving the Poisson equation. 24 
Figure 0] shows that in a PbTiC>3 slab of two unit cells, 
where A2/A1 ~ 1.1, metallic screening from the elec- 
trodes results in a potential that is close to constant in 
the PbTiC>3 slab. On the other hand, in a BaTiC>3 slab 
where A2/A1 ~ 0.9, a similar calculation yields a screen- 
ing potential with a significantly larger depolarizing field. 
These Thomas-Fermi results can be explained as follows. 
Suppose that the metal screens a fraction / of the poten- 
tial drop A2. For A2 > Ai, this translates into screening 
a larger fraction /A2/A1 of the ferroelectric potential 
drop Ai. However, if A2 < Ai, then the fraction of 
the ferroelectric potential drop, /A2/A1, is less than /. 
The difference in the polarization stability of ultrathin 
PbTiC>3 and BaTiC>3 films is therefore directly related 
to the relation of A2 and Ai for the two materials. We 
therefore emphasize the importance of the inequivalence 
of the work functions in determining the ferroelectric be- 
havior of ultra-thin films at this range of thickness. 
In this paper, we have presented the first ab initio 



demonstration of realistic electrodes stabilizing polariza- 
tion normal to the surface in ultra-thin films of < 10 A 
thick. We have shown that proper electrical and chemi- 
cal boundary conditions are essential in stabilizing ferro- 
electricity. Using various electrodes, we find that mon- 
odomain fcrroelectricity can persist down to one unit cell, 
suggesting that thin film ferroelectrics are not specific to 
a single system. The ferroelectric polarization results in 
a difference in the work functions at oxide surfaces that 
must be considered in modeling the metallic screening. 
We demonstrate this difference in PbTiC>3 and BaTiC>3 
films in which we find significantly different screening be- 
haviors with the same electrodes. 
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